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Band structure calculations for Ba2Cu304Cl2 within the local density approximation (LDA) are 
presented. The investigated compound is similar to the antiferromagnetic parent compounds of 
cuprate superconductors but contains additional Cus atoms in the planes. Within, the LDA metal- 
lic behavior is found with two bands crossing the Fermi surface (FS). These bands are built mainly 
from Cu 3dj.2_y2 and O 2px,y orbitals, and a corresponding tight-binding (TB) model has been pa- 
rameterized. All orbitals can be subdivided in two sets corresponding to the A- and B-subsystems, 
respectively, the coupling between which is found to be small. To describe the experimentally ob- 
served antiferromagnetic insulating state, we propose an extended Hubbard model with the derived 
TB parameters and local correlation terms characteristic for cuprates. Using the derived parameter 
set we calculate the exchange integrals for the CU3O4 plane. The results are in quite reasonable 
agreement with the experimental values for the isostructural compound Sr2Cu304Cl2. 
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Several new cuprate materials were studied in the last 
years which have modifications in the standard Cu02 
plane. That class includes not only the now fa- 
mous ladder compounds Sr„_iCu„+i02n ffl, but also 
Ba2Cu304Cl2 or the similar Sr2Cu304Cl2 The 
latter two compounds are antiferromagnetic insulators 
which are characterized by a CU3O4 plane with addi- 
tional Cu_B atoms in the standard Cua02 plane (Figs. 1 
and 3). This new class of cuprates is currently of large 
scientific interest since it may shed light on some of the 
open questions in high-Tc superconductivity. The lay- 
ered cuprate Ba2Cu304Cl2 itself is also interesting due 
to its rich magnetic structure and the unconventional 
character of its lowest energy excitations. Experimen- 
tally, two Neel temperatures have been found ~ 330 
K and ~ 31 K connected with the two sub- 
lattices of A- and B-copper. The magnetic susceptibility 
and the small ferromagnetic moment have been explained 
phenomenologically together with a determination of 
the exchange integrals. Like in undoped Sr2Cu02Cl2 d), 
the lowest electron removal states in Ba2Cu304Cl2 can 
be interpreted in terms of Zhang-Rice singlets with 
a new branch of singlet excitations connected with the 
B-sublattice Q . 

One might expect that the additional Cus give rise 
to considerable differences in the electronic structure in 
comparison with the usual Cu02 plane. In particular, 
the amount of coupling between both subsystems seems 
to be crucial. That question shall be addressed here by a 
microscopic investigation starting from a bandstructure 
calculation. However, as for most undoped cuprates, the 
local density approximation (LDA) leads to a metallic 



behavior and one has to treat the electron correlation in 
a more explicit way. In this paper, we use bandstructure 
calculations to determine tight-binding (TB) parameters 
and then add the local Coulomb correlation terms. In 
contrast to the preliminary TB fit |jll|] we now distin- 
guish between different oxygen orbitals. Furthermore, 
we estimate the exchange integrals. 
A. Bandstructure calculation. The tetragonal unit 
cell of Ba2Cu304Cl2 is shown in Fig. 1. We performed 
LDA calculations for this substance using the linear com- 
bination of atomic orbitals (LCAO). Due to the relatively 
open structure 8 empty spheres per unit cell (between 
Cua and CI) have been introduced. The calculation was 
scalar relativistic and we have chosen a minimal basis 
consisting of Cu(4s,4p,3c?), 0(2s,2p), Ba(6s,6p,5d) and 
Cl(3s,3p). To optimize the local basis a contraction po- 
tential has been used at each site [|l2| . The Coulomb part 
of the potential was constructed as a sum of overlapping 
spherical contributions and the exchange and correlation 
part was treated in the atomic sphere approximation. 
The resulting bandstructure is shown in Fig. 2a. Its main 
result are two bands crossing the Fermi surface (FS). 
These two bands, and most other bands also, have nearly 
no dispersion in z-direction which justifies the restriction 
to the CU3O4 plane in the following discussion. Analyz- 
ing the orbital character of the two relevant bands we 
found that the broad band is built up of predominantly 
the Cua 2idx2_y2 orbital which hybridizes with one part 
of the plane oxygen 2pj;^y orbitals resulting in dpa bonds. 
The corresponding oxygen orbitals, directed to the Cua 
atoms, will be denoted here as p orbitals. They are to 
be distinguished from the oxygen tt orbitals which are 
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perpendicular to them (see Fig. 3). Those oxygen tt 
orbitals hybridize with Cub 3dx2-y2 building the narrow 
band at the FS. So we have to consider 11 orbitals in 
the elementary cell of CU3O4, namely 2 Cu^ Sd^^-y^, 1 
Cus 3^2.2 _j,2, 4 oxygen p and 4 oxygen tt orbitals. In 
Fig. 2b we pick out the corresponding bands from the 
LDA-bandstructure for which the sum of all 11 orbital 
projections is large. It is seen that the two bands cross- 
ing the FS have a nearly pure 3dx2-y2 and 2px,y character 
. It is important to note that the sum of all orbital 
projections in Fig. 2b decreases with increasing binding 
energy. That can be explained since for larger binding 
energy more weight goes into the overlap density in our 
LCAO which does not contribute to Fig. 2b. One can 
also observe that the lower 8 bands in Fig. 2b are not as 
pure as the upper three. For the upper bands only a very 
small weight of additional orbitals, in particular Cu 4 s 
contributions, has been detected. These contributions 
are neglected in the following. 

B. Tight-binding parameters. It is our main goal to 
find a TB description of the relevant bands crossing the 
FS. This task is difficult due to the large number of bands 
between -1 eV and -8 eV of Fig. 2a. There is no isolated 
band complex which makes a TB analysis easy. But we 
have pointed out already that the relevant bands are a 
nearly pure combination of Cu dx2-y2 and O Px.y orbitals. 
So we will only concentrate on these orbitals, accepting 
some deviations in the lower band complex between -3 
and -8 eV. The relevant orbitals are depicted in Fig. 3. 
We distinguish two classes of orbitals, one consists of Cu^ 
3dx2^y2 orbitals with on-site energies and oxygen p or- 
bitals with Ep, the other class incorporates Cus 3^3,2 _j,2 
with on-site energies and oxygen tt orbitals with e^r. 
The coupling between both classes of orbitals which cor- 
respond to the Cvla- and Cus-subsytems, respectively, 
is provided by the parameter tp-jr- Corresponding to the 
different orbitals, one has to distinguish 4 on-site energies 
(e^, e^, Ep, Ett), the nearest neighbor transfer integrals 
(tpd, ^Trd), and several kinds of oxygen transfers (tpp, t\^, 
i^^, ipTr), all together 10 parameters, which are sketched 
in Fig. 3. 

We found that the Cua-Cub transfer tdd can be neglected 
since its estimation yields a value smaller than 0.08 eV 
. Each p orbital is located between two Cua sites and 
we neglect the influence of different local environments on 
the tpp transfer integral. In the case of the ^tttt transfer 
we distinguish the two possible local arrangements, but 
the numerical difference between and t^^ is small (see 
Table ) . The necessity to distinguish between oxygen p- 
and TT-orbitals was first pointed out by Mattheiss and 
Hamann ||l3| for the case of the standard Cu02 plane. 
Since there is a considerable admixture of other orbitals, 
especially Cu idxy and Cu 'id^z^_r2 , in some of the lower 
bands of Fig. 2b, we cannot determine the 10 TB parame- 
ters by a least square fit of the 11 TB bands to the heavily 
shaded LDA bands of Fig. 2b. Instead, at the high sym- 



metry points r = (0,0) and M = {n/a^w/a) we picked 
out those bands in Fig. 2b which have the most pure 
3dx2-y2 and 2px^y character. Only those energies were 
compared with the TB bandstructure (Fig. 2c) derived 
by diagonalizing a 11 x 11 matrix. In this way it is pos- 
sible to calculate the parameter set analytically because 
the TB matrix splits up into 3x3 and 4x4 matrices at 
the high symmetry points T = (0, 0) and M — (tt/o, n/a). 
The resulting parameters are given in Table . These val- 
ues are similar to that which are known for the standard 
Cu02 plane. The largest transfer integrals are tpd = 1.43 
eV and t^^d — 1-19 eV as expected. But they are some- 
what smaller than in the previous TB fit where all 
oxygen orbitals have been treated to be identical. The 
difference between tpp and i^rTr is roughly a factor of 2 
in coincidence with the situation in the standard Cu02 
plane |]l3|. We have found that only the smallest pa- 
rameter, tpT^r = 0.25 cV, is responsible for the coupling 
between the subsystems of Cua and Cub. Thus despite 
the fact that the two oxygen p and tt orbitals are located 
in real space at the same atom, they are quite far away 
from each other in the Hilbert space. 
C. Exchange integrals. Thus far we have found that 
the TB parameters are rather similar to the standard 
Cu02 case and that the coupling between Cu^- and 
Cus-subsystem is quite small. This justifies the usage 
of standard parameters for the Coulomb interaction part 
of the Hamiltonian. Of course it would be desirable to de- 
termine these values by a constrained density functional 
calculation for Ba2Cu304Cl2, but we expect only small 
changes in the estimation of exchange integrals presented 
below. 

The Coulomb interaction also changes the on-site copper 
and oxygen energies. Their difference, given in Table , is 
too small to explain the charge transfer gap of ~ 2 eV in 
Ba2Cu304Cl2 ll6|. Adding 2 eV to the on-site oxygen en- 
ergies, the difference A = Ep — (in hole representation 
which is chosen from now on) becomes similar to the stan- 
dard value derived by Hybertsen et al. jl^ for La2Cu04. 
We have used the values of Ref. [0 also for Ud, Up, Upd 
and Kpd- Since we now have two oxygen orbitals at one 
site we also have to take into account the corresponding 
Hund's rule coupling energy. Unfortunately, that corre- 
lation energy is known with less accuracy than the other 
one and we choose here the simple rule Jj^ — — O.lC/p. 
The Coulomb repulsion between two oxygen holes in p- 
and TT-orbitals is assumed to be Upi^ = Up + 2Jj^, which 
is a valid approximation given degenerate orbitals. In 
the second line of Table we combine the TB parameters 
derived from the bandstructure of Ba2Cu304Cl2 (now in 
hole representation) with the standard Coulomb correla- 
tion terms. This parameter set then defines a 11 band 
extended Hubbard model for the CU3O4 plane which is 
used for the following estimation. 

The exchange integrals have been calculated using the 
usual Rayleigh-Schrodinger perturbation theory on small 
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clusters (Fig. 4). All transfer integrals (and jfj^) have 
been considered as a perturbation around the local limit. 
We calculated all exchange integrals in the corresponding 
lowest order. The exchange Jaa oc t^^/A^ between two 
Cu^ spins is given in 4th order for the simple Cu^i-O- 
Cu-A cluster. It turns out that the influence of intersite 
Coulomb and exchange terms Upd and Kpd is rather large, 
decreasing Jaa from 246 meV to 99 meV (see Table ). 
In spite of our rather approximate procedure, the latter 
value agrees quite reasonably with the phenomenological 
value (130 ± 40) meV § for Sr2Cu304Cl2. Jaa is thus 
also quite close to the standard value of the Cu02 plane 
(~ 140 meV 0). 

The exchange Jb b is given only in 6th order for a larger 
cluster of two Cub, one Cvla and 4 oxygen orbitals. 
Correspondingly, it is roughly one order of magnitude 
smaller, Jbb ~ 17 meV (Table ). For Jab we need to 
distinguish antiferromagnetic and ferromagnetic contri- 
butions. There are two AFM couplings between nearest 
neighbor copper atoms J^^^ and third nearest neigh- 

bor copper atoms J\g ^ j , both being comparably small 
at 4.6 and 0.8 meV, respectively. The ferromagnetic con- 
tribution Jabj— -20meV between nearest neighbor cop- 
per spins arises in 5th order and is provided by Hund's 
rule coupling of two virtual oxygen holes sitting at the 
same oxygen. Due to the uncertainty in this value 
has to be taken with care. 

Summarizing, we presented a LCAO-LDA bandstructure 
calculation for Ba2Cu3 04Cl2. Deriving TB parameters 
from it, we found only a weak coupling between the two 
sets of orbitals connected with the subsystems of copper 
A (Cu^ d^2_y2 and O p) and copper B (Cus d^2_y2 and 
O tt). Furthermore, we found exchange integrals Jaa^ 
Jbb and Jab in reasonable agreement jl^ with phe- 
nomenologically derived values from magnetic suscepti- 
bility data if we add to the TB parameters the standard 
local Coulomb correlation energies. 
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with Upd, Kpd 
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7(3) 
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6.9 
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99 
17 
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-20 
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10 ± 1 

-12 ± 9 



TABLE IL 
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FIG. 1 The body centered tetragonal unit cell of 
Ba2Cu304Cl2 with lattice constants a=5.51 A and 
c=13.82 A . 

FIG. 2 a) LCAO-LDA band structure in the CU3O4 plane 
of Ba2Cu304Cl2, the Fermi level is at zero energy, b) 
The same as in a), but the weight of the lines is scaled 
with the sum of all 11 orbital projections that are used in 
the TB model, c) The band structure of the TB model. 
The parameter set used is shown in Table . The wavevec- 
tor is measured in units of (7r/a,7r/a). 

FIG. 3 Two elementary cells of the CU3O4 plane in 
Ba2Cu304Cl2 or Sr2Cu304Cl2 and the Cu 3dx2-y2 and 
O 2px,y orbitals comprising the TB model. Also shown 
are the corresponding transfer integrals tpd, tTrd, tpp, t\^, 
t^^ and . The Cua orbitals with onsite energy are 
marked by black diamonds, the Cus orbitals with onsite 
energy by black squares and the two different kinds of 
O orbitals with onsite energies £p and e^, respectively, by 
black circles. The orbitals of the B-subsystem are shaded 
to distinguish them from the orbitals of the A-subsystem 
(white) . 

FIG. 3 Clusters used for the calculation of the exchange 
integrals a) Jaa, b) Jbb, c) Jab^/ ^^^d J^bj, d) 

(3) 

"^AB af- Cvla sites are marked with A, the Cub 

sites with B and the oxygen sites with O. 

TABLE I Parameters of the TB fit and the pro- 
posed extended Hubbard model for the CU3O4 plane in 
Ba2Cu304Cl2. 

TABLE II Different exchange integrals as explained in 
the text. Compared are estimations within the extended 
Hubbard model with experimental values m. 
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Figure 1: Rosner et al. 



Figure 3 Rosner et al. 




Figure 4 Rosner et al. 
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